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SUMMARY 


Research carried out under NASA grant NAG 5-145 (Including Its 11 
supplements) has produced a theory which predicts tides in turbulent, 
self-gravitating and loading oceans possessing linearized bottom friction, 
realistic bathymetry and continents; at coastal boundaries "no-flow" 
conditions are imposed. The theory is phrased in terms of spherical 
harmonics, which allows the tide equations to be reduced to linear matrix 
equations. This approach also allows an ocean-wide mass conservation 
constraint to be applied. Solutions have been obtained for 32 long- and 
short-period lunl-solar tidal constituents (and the pole tide), including the 
tidal velocities in addition to the tide height. Calibrating the intensity of 
bottom friction produces reasonable phase lags for all constituents; however, 
tidal amplitudes compare well with those from observation and other theories 
only for long-period constituents. 

In the most recent stage of grant research, traditional theory (Liouville 
equations) for determining the effects of angular momentum exchange on Earth’s 
rotation were extended to encompass high-frequency excitations (such as 
short-period tides). This required incorporating a frequency-dependent 
response of the oceans to the rotational perturbations Induced initially by 
those excitations, as well as frequency-dependent core decoupling. 
Determination of the oceanic responses, and therefore actual calculation of 
short-period tidal effects on Earth’s rotation, would not have been possible 
without the use of the tide theory previously developed by the author. Such 
calculations were carried out for a variety of Initial excitations, Including 
the 32 tide solutions previously obtained by the author and also a number of 
tide models determined by other theorists and from satellite altimetry — with 
tide velocities of the latter models computed by a combination of my theory 
and the models’ tide height harmonics. 

A number of graduate student research projects were also supported by 
grant NAG 5-145. These projects helped to improve our theoretical 
understanding of the dynamic pole tide in the North Sea; the equilibrium pole 
tide world-wide, and its effects on the Chandler wobble; and the effects of 
long-period luni-solar ocean tides on the length of day. 
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IF 0 INAL TECHNICAL REPORT 


The common, underlying theme of all research projects I have carried out 
under support of NASA grant NAG 5-145, including its eleven supplements and 
extending from March 1981 through December 1991, has been the role of the 
oceans in problems of central Importance to geodesy. The primary achievement 
of my grant research — see the frontispiece of this proposal as well as the 
discussions below — was the development of a sophisticated theory for the 
prediction of tides in mid-ocean; its primary applications were to extend 
rotational theory to include excitations of polar motion and UT1 at 
frequencies as high as semi-diurnal, and to accurately predict ocean tidal 
effects on rotation using the extended theory. 

This final technical report for Grant NAG 5-145 will be considerably 
shorter than one might consider appropriate for a final report encompassing 
more than 10 years of funded research. I believe such brevity is acceptable 
because a very detailed and comprehensive report had been submitted two and a 
half years ago, at the request of incoming NASA program managers. 


REVIEW OF RESEARCH ACHIEVEMENTS DURING GRANT PERIOD 


I. My ocean tide theory was initially developed in the context of the 

pole tide, the oceans’ response to the Chandler wobble, in order to quantify 
the contributions of the oceans to the Chandler wobble period and decay rate. 
Originally, however, my grant research began with a more novel approach to the 
oceanic effects on polar motion. The idea was to treat the oceans as an 
integral body, capable of independent rotation except for the torques exerted 
on it by the solid earth which partially contains it. The oceans thus Impart 
an extra "degree of freedom" to the rotating Earth. By analogy with the 
oscillations of a pair of coupled springs or pendulums, the coupled, rotating 
ocean/solid-earth system should possess at least two natural wobble modes. 

With strong coupling (and, of course, the oceans are almost perfectly coupled 
to the solid earth), one wobble mode of the system will correspond to the 
Chandler wobble — yielding incidentally the oceanic contribution to the 
Chandler period. At the time this project was proposed, it was anticipated 
that the frequency of the second mode would be close to the first (for very 
strong coupling), thereby explaining the observed modulation ("beating") of 
the Chandler wobble amplitude. 

The theory was worked out [Dickman 1983] for the case of an elastic 
mantle possessing a fluid core, with fluid oceans that respond statically to 
the incremental centrifugal force incurred from wobble. The torques acting 
between the oceans and solid earth were parameterized in terms of the relative 
angular velocity between the two bodies, and the wobble mode frequencies, 
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decay rates, and relative amplitudes were determined for a wide range of 
coupling strengths. Unexpectedly (though explainable in hindsight), it turned 
out that — because of the shape and fluidity of the oceans — the second 
wobble mode was necessarily retrograde (clockwise wobble motion); at any 
coupling strength it would not lead to a second Chandler frequency, since that 
wobble is prograde. Surprisingly, at high but physically reasonable coupling 
strengths, the second mode exhibited many characteristics of the controversial 
"Markowitz wobble", a low-amplitude ultra-long period apparent polar motion 
[see, e.g., Dickman 1981] which many researchers claim to be pure noise. 

Thus, the results from the first grant project appeared to indicate that an 
Earth with oceans possesses two natural wobble modes rather than Just one. 

The theory developed for that project, however, was based on an assumed 
equilibrium behavior of the fluid oceans; it was argued that this 
approximation would not be far wrong because the wobble periods are so long 
(in excess of a year). Such arguments were not universally accepted; it is 
not, after all, intuitively obvious how a wobble of the oceanic "body" could 
survive fluid dynamic motions [cf. Wahr 1984, Dickman 1985a]. 

The extent of fluid dynamic disruption of ocean/solld-earth wobble at 
physically reasonable coupling strengths is related to the extent to which the 
oceanic response to wobble departs from static; that is, it is related to the 
extent to which the pole tide is non-equilibrium. Observations of the pole 
tide for the past century have, almost without exception, indicated it to be 
significantly enhanced — at least in a number of shallow seas, and probably 
in the deep ocean as well — and capable of significant contributions to the 
Chandler wobble period and decay rate [see Naito 1979, Dickman 1979, Lambeck 
1980]. The expected low amplitude of the pole tide (~ 0.5 cm for typical 
wobble amplitudes) and relatively high noise level of most tide data, however, 
cast doubt on the accuracy of the data and on the statistical reliability of 
even the best analyses. [But in time series analysis, long durations of data 
(and repetitions of the periodicity in question) always allow a low-amplitude 
signal to be picked out. It turns out (see below) that most likely the data 
analyses were not wrong, but misinterpreted: they were revealing coastal 
features of the pole tide, unrelated to its global character.] Clearly, it 
was time for a fluid dynamic theory of the pole tide. 


II. Prior to my beginning the development of a dynamic pole tide theory, 

work by others [see, e.g., Munk & MacDonald 1960, Smith & Dahlen 1981, Carton 
& Wahr 1983 (based on a 1981 paper)] had suggested that dynamic effects were 
likely to be very small. Such conclusions had major Implications for the 
extent of low-frequency mantle anelasticity, for core viscosity, and even for 
the evolution of the lunar orbit; they would also resolve long-standing 
questions about the nature of the Chandler wobble. Those conclusions, 
however, depended on approximations either to the equations or to the methods 
of solution; if the pole tide data analyses were indeed correct, it appeared 
that an exact fluid dynamic theory would be necessary. 

I chose a spherical harmonic approach for two reasons: spherical 
harmonics are appropriate for a global description of tidal quantities; and, 
the effects of the pole tide are primarily the result of oceanic mass 
redistribution, i.e. changes in Earth’s inertia generated by the pole tide, 
and the oceanic products of Inertia could be easily represented as a single 
spherical harmonic [see Dickman 1985b]. Thus, the tide height T and 
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north-south, west-east tide current velocities u , u. were expressed as 
complex variables, periodic in time with frequency <r. 


T = Re{7} 

T = T exp ( i<rt ) 


u 9 = Re % } 

A 

U 0 = u 0 exp(i<rt) 


u^ 




= Re{u^> 

A 

= u^ exp(i<rt) 


the time-independent portions of these variables were expanded in complex 
spherical harmonics, 





where is the harmonic function (of colatitude 0 and longitude X) of degree 
l and order n. The sets of unknown coefficients 


t = <i£> iJ = (u£> = (vj> 

are the quantities our theory must determine. 


The frequency <r of the tide in this case would be that of the Chandler 
wobble, which causes the pole tide. However, the period and decay rate of the 
Chandler wobble are partially the result of pole tide effects (the 14-month 
Chandler period would be ~1 month shorter, and the wobble energy would be 
dissipated more slowly, if the Earth were oceanless). This "feedback" 
situation implies that the pole tide characteristics — which depend on <r — 
cannot be determined Independently of <r itself. For this reason, and because 
no approximations would be made in the development, I began the project with a 
simple ocean model: non-turbulent, global oceans of uniform depth possessing 
bottom friction with constant drag coefficient. 


The tide-governing equations, Including two horizontal "momentum" 
equations and one equation of pointwise mass conservation ("continuity"), will 
be presented in full later in this report. For the relatively simple ocean 
model employed at this stage of work, the momentum equations could be solved 
analytically for the tide velocities in terms of the tide height, then 
substituted into continuity to yield a single differential equation for the 
tide height T. With T expanded in spherical harmonics, it became necessary to 
formulate a general expression for derivatives of harmonic functions; after I 
developed the required relations [see Dickman 1985b for details] it became 
apparent that the theory would depend on the computation of triple harmonic 
product integrals, viz. 


nqa 

tfxi 


J 


(Y?)*yV^ ds 
l fX j. 


(these integrals are over the unit sphere). Eventually it was possible to 
reduce the tide-governing equations to 


(1) 


B • f = M t , 

= “ P * 


where the polar motion forcing the tide is of amplitude M . For the ocean 
model considered at this stage the components of equatlon P (l) are relatively 
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uncomplicated; for example, only the tide height coefficients T^ are non zero 
and need be determined [see Dlckman 1985b for details and properties of B and 


Joint solution of equation (1) and the Liouville equation, the latter 
describing conservation of angular momentum, allowed the pole tide (T) in 
global oceans, and its effects on polar motion (cr), to be determined 
self-consistently. The computed solutions included the effects of both the 
tidal inertia (associated with changes in sea-level) and tidal relative 
momentum (associated with tidal currents). Overall, it was found that the 
pole tide in global oceans is very close to equilibrium; however, the slight 
dynamic effects on period and decay time Implied by these results could well 
be amplified by other factors, such as non-globallty or self-gravitation, in 
real oceans. In short, the ocean model used in this stage of my research was 
simply too preliminary to settle the question of pole tide effects on wobble. 

The tidal theory was extended to include more realistic ocean models in 
later stages of the grant research. Such work is described in the sections 
following the next. 


HI. One region where the pole tide has been consistently observed with 

greatly enhanced amplitudes — up to 10 times equilibrium — is the North and 
Baltic Seas area. Such enhancements, in those seas alone, would be sufficient 
to dissipate up to half of the Chandler wobble energy [Dickman 1979; see also 
Dlckman 1986). Given the apparently nearly static behavior of the pole tide 
globally, it was important that such shallow-sea observations be verified 
theoretically. 

A fluid dynamic model of the North Sea pole tide was indeed proposed in 
the mid-1970’ s [Wunsch 1974], and it appeared able to explain the observations 
— in particular, its combination of depth-dependent bottom friction and 
shallowing to the south generated an "eastward Intensification" of the tide 
height similar to what had been observed along the North Sea coast. However, 
Wunsch [1975] later declared, with little elaboration, that his results were 
erroneous . 

In the earliest years of the grant period, I had supervised a graduate 
student (J. Prelsig) who re-evaluated the analysis of Wunsch in light of his 
later erratum. Our review of Wunsch [1974 and 1975] had revealed that his 
solutions exhibited a mild eastward intensification; but they failed to 
satisfy one of his tidal equations and either did not satisfy fundamental 
boundary conditions or else depended on invalid approximations. We found that 
a slightly different model of the North Sea bathymetry avoided such 
difficulties. Our better-quality solution would later allow me to discover 
the essential problem with North Sea pole tide dynamics. 

I began my grant research for this project by re-deriving the North Sea 
tide equations, paying special attention to the various approximations Invoked 
in earlier works. Obtaining (or searching for) solutions to these equations 
required an extensive effort. All solutions [see Dickman & Preisig 1986] 
shared several noteworthy features: 1) they all exhibited an eastward 
intensification along the south shore of the Sea — more so than Wunsch’ s 
solutions, and in good agreement with observation; 2) they all dissipated far 
more Chandler wobble energy than is physically possible (with the most Intense 
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currents occurring In the northern areas); and 3) they all were non-unique 
solutions. The Implication was that the North Sea tide equations themselves 
were deficient. 

The non-uniqueness and extreme Intensity of the solutions bear a 
resemblance to the elements of Stokes' Paradox’, details may be found In 
Dlckman & Preislg [1986]. Like Stokes’ Paradox, these features suggested that 
the inclusion of inertial terms in the North Sea tide equations would produce 
a more reasonable solution. The Inertial terms would either be time 
derivatives (3u/3t for the momentum equations, 9T/flt for continuity), which 
are normally a part of "Laplace" tide equations but were omitted by Wunsch 
[1974]; or — more likely — advectlve terms (u’Vu In the momentum equations), 
which are normally excluded from tide equations because they are non-linear, 
but which might be important in shallow, partially confined seas. 

A rigorous fluid dynamic theory of the pole tide, including inertial 
terms, has yet to be worked out for the North Sea; thus, the observations of 
markedly amplified pole tide heights In the North and Baltic Seas have yet to 
be theoretically substantiated. In light of the work described in the next 
section — which demonstrates that the dissipation of wobble energy by the 
deep-ocean pole tide is negligible — the North-Baltic Seas observations alone 
stand in the way of attributing the bulk of Chandler wobble dissipation to 
mantle anelastlcity. Without solution of this problem, our knowledge of 
mantle rheology at low frequencies remains speculative. 


During this time I had Just finished supervising another graduate student 
(D. Steinberg). He had (among other things) analyzed tide data at a number of 
ports — mid— ocean, coastal, and shallow-sea; comparing the results with those 
from carefully constructed equilibrium pole tide time series, and accounting 
for the high noise level in the tide data, he showed that the pole tide in 
mid-Pacific (Hawaii) and mid-Atlantic (Bermuda) was statistically 
significantly greater than static, perhaps 1.5 to 2 times greater. The 
calculations of the static pole tide time series, incidentally, were based on 
an extension I developed of earlier theory using suggestions by Dahlen [1976]; 
that development led me to appreciate the distinction between "untruncated" 

(or global ) and "truncated" (i.e. restricted to the oceans ) tidal quantities 
[see Dlckman & Steinberg 1986 for details]. That distinction would play a 
major role in my subsequent work. 

The potential effects on the Chandler wobble of a global pole tide with 
the characteristics suggested by our careful analysis — see Dlckman & 
Steinberg [1986] — were severe enough, and their consequences for mantle and 
core properties were important enough, that further theoretical work on pole 
tide dynamics in realistic oceans was warranted. 


IV. The importance of accounting for sea-floor topography and land 

barriers (i.e. continents) in tidal theory is matched only by its difficulty. 
The extension of my global dynamic pole tide theory to the case of non-global 
oceans faced several obstacles. First, a variable ocean depth h makes the 
tide equations more complicated; with h = h(0,A) expanded in spherical 
harmonics (h - £ h. Y^), for example, the development would Involve quadruple 
rather than triple^pr^ducts of harmonic functions. Second, an adequate 
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spherical harmonic description of oceanic bathymetry had been unavailable, and 
would need to be constructed. Third, the existence of oceanic boundaries 
requires that the tidal currents satisfy no-flow conditions at the boundaries. 
Fourth, the non-globality of the oceans (i.e. the presence of continents) 
changes the character of the polar motion, imparting an ellipticity to the 
(otherwise circular) Chandler pole path. This last difficulty implied that 
the polar motion, of complex amplitude m = m + 1m , will contain prograde and 
retrograde components: X ^ 


( 2 ) 


m 


i<r t 

= M e - 


M e 

~R 


-i<r*t 


Since the potential which forces the pole tide depends on m, equation (2) 
implies that the tidal potential will be more complicated as well. 


These problems were treated as follows [see Dickman 1988a for details] . 


Y ±1 Y^, 


1. Relations which allow products of harmonic functions (e.g.Y” Y . , 
and so on) to be reduced to sums of harmonic functions were ' 

laboriously developed, then selectively applied to reduce quadruple products 
to triple products, for which the integrals (A^J^, see earlier) are easily 
computed. It was then possible to re-write the^tide-governing equation as a 
single non-differential matrix equation of the form 


(3) 


B • ? - M t> + M*S 

= ~P P “R R 


where now the collection of unknown tide height coefficients f = {T^} involves 
all degree and order harmonics. The elements of the matrix B and vectors b p , 

, which are significantly more complicated than in the case of global oceans 
(--it now took ~ 5 sheets of notepaper to write them out symbolically), depend 
on the tidal frequency bottom friction P, and the oceanic bathymetry 
coefficients h.. 


2. At the time this research was performed, published values of h. 
were available only through degree and order 8. Given the other obstacles to^ 
this research, I decided to postpone construction of a more complete set of 
bathymetry coefficients (see next section). From the work Reading to Dickman 
& Steinberg [1986], we had obtained two extensive sets of 0., the coefficients 
of the ocean function (see below), so that an accurate description of the 
continent - ocean distribution was already in hand. Thus, the results at this 
stage of my work would represent tides in non-global but flat-bottom (i.e. 
uniform depth) oceans. 


3. The constraint that continents are barriers to tidal currents — 
that is, there is no flow of tidewater across continental boundaries, is 

u • n = 0 at coastlines 


where n is the unit normal to the coast and u = (u^,u. ) is the horizontal tide 
velocity. I re-wrote this "pointwise" constraint as the more suitable, global 
constraint 

(4) ill • ^(h Q 0) = 0 everywhere 

where h is a constant (for scaling the boundary conditions) and 0 is the 
ocean function, defined as being zero at land locations and unity over the 
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oceans. With 0 expanded In spherical harmonics, and u^, expressed In terms 
of the tide height T, It was possible to reduce (4) to” a single 
non-dlfferentlal matrix equation, 


(5) 


IB' • f = M t' + . 

= "P P HR 


The pole tide in non-global oceans, constrained at coastlines, is 
determined as the joint solution to equations (3) and (5). Eventually, the 
approach I adopted was to use an equal number of tide equations (comprising 
(3)) and boundary condition equations (comprising (5)), thus altogether a 
greater number of equations than unknowns: an M overdetermined" situation; the 
actual matrix inversion procedure to obtain f would yield a least-squares 
solution — given the factor h Q , a weighted least-squares solution. 

4. Non-global oceans modify both the wobble frequency (a) and the 
proportion of retrograde to prograde wobble motion, Z * M*/M . I found a way 
to rewrite the Liouville equations in this case so that <r K and Z could be 
determined, once M had been prescribed and a tide solution obtained; see 
Dickman [1988a] for details. 


Thus, the self-consistent dynamic pole tide in non-global oceans would be 
found as the joint solution to Liouville equations and equations (3) & (5). 

The complexity of the problem, and need to search for joint solutions, meant 
that a supercomputer (like the CYBER 205 or IBM 3090 vector facility) was 
required. Early results suggested that non-globality had a tremendous effect 
on the Chandler wobble [see Dickman 1988a for details]. However, a number of 
theoretical tests eventually led me to realize that my solutions generate a 
non-zero tide on land, as an automatic artifact of analytical tide equations; 
that only the oceanic portion OT of the solution should be employed (e.g. in 
the Liouville equations) to determine tidal effects on rotation; and that the 
land portions of the solution could safely be ignored. 

The “truncated” harmonics of the solution, i.e. the coefficients of OT, 
were easily computed from the solution for f. The results indicate that pole 
tide dynamical effects on wobble are small. For example, tide dynamics add 
5/3 days more to the Chandler period than a static pole tide would, but 
dissipate a negligibly small amount of wobble energy. 

During this stage of research I also modified the non-global theory to 
account for the response of the elastic mantle to the tidal force, and to 
account (at this point somewhat approximately) for self-gravitation and 
loading by the tide layer. Table 3 in Dickman [1988a] (not presented here) 
summarizes the pole tide effects on wobble including these factors, for a 
range of bottom friction strengths. Although loading and self -gravitation 
increase the height of the ocean tide, the corresponding depression of the 
sea-floor reduces the solid earth's contribution to the Chandler period. The 
net effect, as shown in that table, is even less dissipation of wobble energy 
by the pole tide, and a lengthening of the Chandler period by Just over 1 day 
more than that from a static tide. 

At this point in my research, my tide theory was based on an ocean model 
which lacked turbulence, realistic bathymetric variations, and exact loading/ 
self-gravitation. The long period of the Chandler wobble implied that none of 
these deficiencies would be likely to change our conclusions significantly 
[cf. Dickman 1989]. For the study of luni-solar ocean tides (almost all of 
which are much shorter period than the pole tide), these factors would of 
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course need to be included for accurate modeling. 


The question of how the oceans respond to atmospheric pressure variations 
has been a recurrent one in wobble and l.o.d. research. With the recently 
recognized importance of atmospheric angular momentum in Earth rotation [e.g., 
Eubanks et al. 1985, Dickey & Eubanks 1986], an answer to that question 
becomes even more crucial: essentially, the extent to which the sea surface 
depresses under barometric high pressures ( — the ideal static response, the 
"Inverted barometer", represents a perfect "lsostatic" compensation) 
determines how much atmospheric pressure excitation of wobble or l.o.d. is 
canceled out. 

My tidal theory could be used to resolve this question. An atmospheric 
pressure fluctuation of amplitude p Is equivalent to a tidal potential 
perturbation of amplitude -p /p [Defant 1961]. However, my theory had to be 
generalized because the barometric forcing could possess any spatial structure 
(and, for that matter, any time-dependence), not Just the degree-2 spherical 
harmonic structure of a tidal forcing. Thus, we write 

P A = Re{p A > B A = (Za^Y^Je*-* 

to describe the atmospheric forcing at frequency o\ 

Letting a = {a?) denote the "vector" atmospheric forcing, I found that 
the equations governing the "tide height" or oceanic response, T = {T^}, could 
be written 

(6) B • f = - — 1_ C B • a . 

P w 8 - 

The matrix B would have been the same as in equation (3), except that at this 
point I had'taken the opportunity to account for loading and self-gravitation 
exactly — a generalization of my earlier tidal theory (atmospheric loading & 
self-gravltatlon were also Included now). The "no-flow" oceanic boundary 
conditions now took the form 


(7) B' • f = °B' • a 

p w* = 

with similar comments about g' In relation to equation (5). See Dlckman 
[1988b] for details regarding B, °B, B' , and °B' . 

Solutions to (6) and (7) were computed for atmospheric forcing described 
by a single harmonic of degree l, order n at period 2n/<r. To evaluate the 
range of oceanic responses, I considered 2n/<r = 5 days, 50 days, and 433.2 
days; for each frequency I tried l=2, h=0; l=2, n=1 ; L=5 , h=5; and L=10, N=10 
(note that l=2 corresponds to global, l=5 to "basin" -scale, and l=10 to 
regional forcing). For flat ocean basins I found that the actual, dynamic 
response does not differ from the Ideal, Inverted barometer response by more 
than 10% - 15%, especially at forcing periods exceeding a week (my work did 
reveal a slight, 9-day resonance, however). The effects of bathymetry and 
turbulence, both of which could be very influential at moderate or short 
periods, were not considered at this stage. 
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V. The addition of turbulence to my ocean model required extensive 

re-worklng of the spherical harmonic theory. The mathematical complexity of 
the turbulent drag forces prevented the three starting equations (2 momentum, 

1 continuity) from being combined Into a single equation reducible to matrix 
form; Instead, each equation had to be written directly as a matrix equation. 
This in turn resulted in the appearance of modified harmonic product 
Integrals, 

J Y^sine ds 

which are similar to but which had never been solved before. Once I had 

worked these out, It was^osslble to write the momentum equations symbolically 
as 

_ZZ7 • 3 + c 3 * ^ 

777 • 3 + ^ ^ . 


relating the unknown tide velocities to the tide height. With these solved 
for 6 and ^ in terms of t, they could be substituted Into the continuity 
equation, yielding a matrix equation for f analogous to equation (3); when 
substituted into the boundary constraint, (4), they would yield a second 
matrix equation for f, analogous to (5). Details of all this may be Inferred 
from Dlckman [1989]. 


The preceding relation between tide velocities and tide height is, to my 
knowledge, the first of Its kind to be developed. When applied to lunl-solar 
tides, for example, it would allow tidal velocities to be determined from the 
various tide height models published by other researchers (whose theories 
yield T but not u , u^). It would also allow tidal velocities to be Inferred 
from SLR or satellite altimetry observations of tide heights. 


The development of the turbulent theory also required that a modified 
ocean depth ^ 

h = h sine 

be expanded in spherical harmonics. At this point in my research I computed 
the values of hu from those of h^ theoretically, for the case of flat ocean 


i 


The results for the pole tide In turbulent, 


h -0(e,X) ). 

self-gravitating & loading oceans possessing bottom friction but 


basins (where h 
non-global , 

no bathymetry were — not surprisingly — very similar to the non-turbulent 
results (except when I employed unrealistically large turbulent viscosity 
coefficients) . 


Around this time I began to construct my own set of spherical harmonic 
coefficients for 0, h, and h directly from actual bathymetry data. The data 
was kindly provided by J. Marsh at NASA/GSFC, In the form of the DBDB5 5’ X 5’ 
digital bathymetry tape. 


11 



VI. The preceding research had answered questions about pole tide 

dynamics in fairly realistic oceans, and about pole tide effects on wobble; 
the work also had implications for dissipation of wobble energy by mantle 
anelastlcity and core viscosity. At this point I was in possession of a 
pretty good pole tide model, and I set about applying it to the case of 
lunl-solar ocean tides. 

The luni-solar tidal potential is of degree 2 order s (|s| =0, 1, or 2 
for long-period, diurnal, or semi-diurnal tides), and has frequency r. We 
wrote it as 

(8) U = Re{ U> , U = MY^ e J - t 


where is the complex forcing amplitude. 

The momentum equations appropriate for the tide in self-gravitating and 
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A, the horizontal eddy viscosity coefficient, governs the strength of the 
turbulent drag forces. These equations account for the effects of Earth 
rotation (f is the Coriolis parameter), the yielding of the solid earth to the 
tide-generating potential (a = 1+k -h is the usual Love-number combination), 
and tidal loading & self-gravitation; the matrix L is defined as 


( 10 ) 
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'lp 


= 5.6 


3 p w l+k'.-h'. 
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21+1 


where k^ , are load Love numbers of degree l, and p is mean Earth density. 

Following the same procedures worked out for the turbulent pole tide, I 
reduced the momentum equations to matrix form and solved for the unknown tide 
velocity coefficients: 
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where E, F, and 0 are defined in Dickman [1989]. As mentioned earlier, in 
situations where tide heights only have been determined, such as satellite 
ranging or altimetric observations, or finite-difference tide theories, 
equations (11) make it possible to estimate the associated tidal currents. 
Such estimates were unobtainable prior to this work. 


( 11 ) 


V 


= E- L- 


-» 

U 


= F- \L- 


The dynamic description of the tides is completed by specifying 
continuity, l.e. polntwise conservation of mass, 

(,2) { slne *[ Ku ° ] • 

With (12) reduced to matrix form and the expressions for the tide velocities 
substituted in, I obtained a single, non-differential matrix equation for the 
tide height, 

(13) B ■ T = Mb . 

— — o 


When the velocity expressions were substituted into a matrix version of the 
global boundary constraint, (4), the result was 





(14) 

h B' • T = h 

M b' . 

0 - 0 

- 0 




b, |', 

and b' , which depend on or, P, A, and 

h^ , are detailed in Dickman 

[1989] . 



With the coefficients of 0, h, and h for the actual oceans constructed 
directly, the solutions to (13) and (14) describe the lunar and solar tides in 
realistic oceans. In Dickman [1989] I obtained such solutions £or 5 long- 
period tides. Comparison of the amplitudes and phases of the T coefficient 
— the primary coefficient of all long-period tides — from my theory with 
those from other theories and observations revealed an overall agreement of my 
amplitudes with the others (but see the discussion in Dickman [1989] for some 
interesting exceptions); however, the phase lags predicted by my theory were 
consistently lower than those found by others; this is discussed later. At 
this time, I also predicted the effects of my long-period ocean tides on 
l.o. d. and wobble; in all cases the effects were found to be small — but 
marginally detectable in recent (and, of course, future) high-quality 
space-geodetic rotation data. 
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VII. One of the major limitations of my o lunl-solar tldaj. predictions was 

an ambiguity resulting from the failure of T to vanish. T relates to the 
total amount of tide water world-wide, p /Ids, and should be zero since tidal 
forces only re-distribute the oceanic mass (but do not add to it). 
Calculations of T and OT are ambiguous depending on whether T° is included in 
their spherical harmonic expansions or excluded (In Dlckman [1989] I had 
found the preferred solutions to be those which excluded T^). 

In fact, other tide theories [e.g. Schwiderskl 1983] deliberately 
generate non-zero T , by allowing their oceans to periodically "leak" water, 
in order to model shallow-sea dissipative processes (Parke [1982] 
additionally arranges the leakiness so that the world-wide net leak is zero). 
A glance at the tables of results in Dlckman [1989] demonstrates how serious 
the ambiguity can be: with the T produced by my theory, for example, tidal 
effects on wobble and l.o.d. gre uncertain by up to 100% in amplitude and 
phase, depending on how the T q is treated. The role of T° was not evaluated 
by other tide researchers, perhaps because their models are not spherical 
harmonic; their analyses thus failed to recognize such ambiguity. 


As discussed in Dlckman [1990], I considered a number of alternatives for 
achieving a reduced (or zero) T°. I Judged the optimal approach to be the 
addition of an oceanic mass-conservation constraint 

f T ds = 0 

ocean 

or, equivalently, 

[CT]° = 0, 

Implemented via Lagrange multipliers. Thus, my best estimate of the tide 
height in realistic oceans was the least-squares solution of 


(15) 
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and A [0T]° = 0 

l o 


where is an undetermined Lagrange multiplier. 


Such mass-constrained tide solutions are the first of their kind. These 
solutions are more realistic, and less ambiguous, than other luni-solar tide 
models. For the five long-period tides I had studied, my predicted amplitudes 
differed noticeably with the constraint. After adjusting the boundary 
condition weighting factor (h Q ) according to other criteria [see Dlckman 
1990], it turned out that the°tide amplitudes were very similar to the 
preferred values in Dlckman [1989] (and thus still in general agreement with 
other estimates). Additionally, the phases showed a slight improvement. 


During this time I also supervised a student (Y.-S. Nam) whose master’s 
project involved tidal effects on the length of day. He constructed a 
relatively long span of high-quality l.o.d. data using the IRIS program’s UT1 
observations; corrected it for atmospheric wind and pressure effects; and 
estimated the 9-day, fortnightly, and monthly tide effects on l.o.d. (the 
analysis was actually phrased in terms of the "zonal response coefficient" 
[Agnew & Farrell 1978]). Using my luni-solar tide program, he computed 
theoretical ocean tide effects at those periods and compared the theoretical 
predictions with observation. 
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His work was carefully researched and he accounted fully for an expected 
core-mantle decoupling. Any statistically significant discrepancies between 
theory and observation would be attributed to factors not Included In our 
model, such as mantle anelasticity. We hoped that, with the differing oceanic 
effects subtracted out, the frequency dependence of those factors could be 
elucidated. Because the effects are small, comparable to the noise level 
present In that VLBI data, it was difficult to draw strong conclusions; our 
results, which are described In Nam & Dlckman [1990], suggest that the 
anelastic contribution to the fortnightly tidal effect on l.o.d. cannot exceed 
1.754, while partial coupling between core & mantle at monthly periods cannot 
be ruled out. 


VIII. It was clear even from my pole tide Investigation [see Dlckman 

1988a] that the tidal response of the oceans at various frequencies depends on 
how close those frequencies are to the natural resonances of the ocean basins. 
For lunl-solar tides, which span three very different frequency bands 
(semi-diurnal, diurnal, and long-period), the dependence on oceanic resonances 
should be even more significant [see Platzman 1984]. The oceanic response to 
tidal forcing must unavoidably be considered strongly frequency-dependent. 

Knowledge of that frequency dependence would be potentially Important for 
the determination of satellite orbits and, simultaneously, satellite 
determination of ocean tides. For example, prior to 1990, the most 
comprehensive estimation of ocean tide constituents from satellite geodesy 
yielded 616 harmonic coefficients of 32 major and minor tides [Chrlstodoulidls 
et al. 1988]; however, the bulk of those coefficients — corresponding to 
minor or "sideband" tides — were merely interpolated from the solutions for 
the major tide coefficients. Chrlstodoulidls et al. [1988] based that 
interpolation on a linear model (l.e. a linear frequency dependence) of tidal 
admittance (tidal admittances are the ratios of actual to static tide height 
amplitude coefficients); although such a model represents an improvement over 
earlier models of constant admittance (In which a minor tide amplitude was 
taken to be In the same proportion to Its forcing potential as the major tide 
In that band), It may not account very well for the effects of resonances In 
or near the tidal bands. A better model of the frequency dependence within 
the tidal bands would improve both the satellite orbit determination and the 
ocean tide Inferences. 

In an effort to elucidate the frequency dependence of tidal admittance, 
and consequently gain insight into the character of natural oceanic 
resonances, I computed the actual oceanic responses to a pre-set forcing, at 
the frequencies of 32 major and minor lunl-solar tides. In addition to 
providing characteristics of tidal constituents (in many cases, unavailable 
prior to this work) of immediate use In satellite geodesy, these responses 
could be used to assess sideband tide models previously employed, and also to 
Infer the existence of oceanic elgenfrequencles. 

In early stages of this research, I took the opportunity to Improve the 
resolution of my spherical harmonic tide model. Using the original (DBDB5 
tape) world-wide togography data, I constructed spherical harmonic 
coefficients of h, h, and 0 through degree and order 48. This took some 
effort, because the computations involved numbers larger than mainframes can 
store, and because the computations had to be performed using extended 
precision. Subsequent numerical tests (and a further look at equation (13)) 


15 



led me to the realization that my programs do not recognize hlgh-resolutlon 
bathymetry unless hlgh-degree tide height coefficients arg al|o obtained 
That Is, the theory automatically and correctly Ignores h. , h. , and ^ 


,n 


harmonic degree j. s 2(l+1) If I solve only for T'^ of 


degree 


0 . of 
L . ^Thus , 


to assess the effects of high- or even moderate-resolution bathymetry I would 
have to run larger programs (on the IBM 3090 vector facility I estimated that 
2-minute, 8-Megabyte programs would have to give way to 35-mlnute, 12-Megabyte 
programs or 4-hour, 32-Megabyte programs in order to Include bathymetry 
exceeding degree 12 or degree 24, respectively). 


Following these modifications I discovered that, with the Increased 
resolution, the spherical harmonic solutions had become sensitive to the 
Intensity of bottom friction (and, marginally, to the extent of lateral eddy 
dissipation). Further experiments revealed that, with a magnitude of P = 
2.5x10" sec" for long-period tides and 1.5x10” sec” for diurnal and 
semi-diurnal tides. It was possible to produce phase lags for the degree 2 
tide height components which were consistent with tide observations and other 
tide theories. This represented a major Improvement In my tide models. 
However, calibrating the Intensity of friction In my model ocean did not 
significantly Improve the abnormally low amplitudes predicted for the 
short-period constituents. 

The results of the study are presented In Dlckman [1991]. One major 
conclusion is that the linear sideband models used in the most recent 
satellite geodesy analyses are definitely better than constant-admittance 
models, but typically still fall to account for ~ 40% of the sideband 
admittance variations. This in turn suggests that theoretically constrained 
satellite analyses might yield much-improved ocean tide models and satellite 
orbit tracks. 


The tidal solutions In Dlckman [1991] were also evaluated with regard to 
a number of "diagnostics" Indicating anomalous oceanic behavior. From those 
clues we concluded that a search for oceanic resonances should focus on the 01 
- Q1 and 001 - J1 frequencies in the diurnal tide band; the M2 - 2N2 
frequencies, and frequencies beyond 2N2, In the semi-diurnal tide band; and 
the MSm - Mm frequencies, and frequencies beyond M9, In the long-period tide 
band. See Dlckman [1991] for details. 


At this time, I also explored the use of external constraints — for 
example, from satellite observations of ocean tides, or from finite-difference 
predictions of tide coefficients — in my spherical harmonic ocean tide 
theory. Initially, it was hoped that, like the "hydrodynamical Interpolation" 
technique used In numerical tide theories, the incorporation of actual data 
would lead to improved short-period tide amplitudes; and, because spherical 
harmonic constraints are global, the results should be less affected by 
coastal anomalies than are the finite-difference theories. 


The constraints, implemented via Lagrange multipliers [cf. Menke 1984] 
according to 


(16) 


f B 

h oS' 


>T = M 
~0 


hfr 


and A [OT] = 0 
l o 


and 


A [OT) = SPECIFIED 
2 L 


where A t and are the undetermined (Lagrange) multipliers for global mass 
conservation and the specified harmonic constraint, respectively, quickly 
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revealed another use: they could be employed to evaluate the fluid 
dynamic validity of tide observations or theoretical predictions. For 
instance, satellite determination of the principal lunar monthly tide height 
coefficient ( [0T]°) by Christodoulldis et al. [1988] yielded an amplitude of 
0.3610.50 cm and" phase lead of 4.1°±76.6°. The surprisingly small amplitude 
(static [0T]°is 1.17 cm, and at such a long period the tide should not be so 
far below equilibrium) as well as the large uncertainties and the phase lead 
suggested that this determination was questionable. To evaluate the satellite 
result, solutions were obtained with [OT] ° constrained to assume the range of 
reported values, also more reasonable values. As described in Dlckman [1991], 
the experiments did Indeed confirm that the satellite value was not realistic. 

Another example considered was the prediction by Schwlderski of the q 
primary semi-annual tide height coefficient, [OT] = (1.24 cm amplitude, 48.3 
phase lag); the excessively large amplitude (~ 20% bigger than static) and 
phase lag of this tide, considering its long period, suggest that his 
prediction is Inaccurate — perhaps contaminated by seasonal meteorological 
effects on the incorporated data. Diagnostics were determined for the tidal 
solution constrained according to Schwlderski’ s value and constrained by 
alternative values. Again, the investigation appeared to confirm that the 
coefficient was fluid dynamically questionable. 

Finally, by using simultaneous multiple constraints it was possible to 
test whether the high- and low-degree harmonics of some tide model are 
mutually consistent. As detailed in Dickman [1991], short-period degree-4 
harmonics as large as were reported by satellite and numerical tide models — 
those admittances were an order of magnitude larger than those in my tide 
model — may well be dynamically incompatible with their degree-2 admittances. 


During the year of that grant project, a graduate student (Y.-S. Nam) 
began to research the problem of wobble of the coupled ocean/solid-earth 
system [Dlckman 1983 — the focus of the first year of this grant]. He was to 
generalize my original work to allow for a dynamic response of the oceans 
during wobble; the dynamic response would be determined using my tide 
programs. Under my supervision, he finally mastered the theory enough to 
begin generalizing it; and completed a first set of computations. The results 
were encouraging but needed checking; however, after writing a summary of his 
efforts, he left the university to undertake doctoral studies elsewhere. 


IX. The goal of my final project under grant NAG 5-145 was the 

determination of tidal effects on Earth’s rotation and on the lunar orbit. 

The tides in question were the 32 long- and short-period constituents I had 
previously obtained as the solutions to my spherical harmonic ocean tide 
model. The effects on Earth’s rotation derive from the change In the oceans' 
inertia tensor, i.e. the redistribution of oceanic mass associated with the 
rise and fall of the sea surface, and from the relative angular momentum 
produced by the tidal currents (see later discussion for details). The 
effects appeared to be calculable from standard Liouvllle equations ( — but, 
again, see later) and after a brief effort I had produced a detailed list of 
the effects of tidal Inertia and relative momentum on both polar motion and 
the length of day. 
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To determine the effects of ocean tides on the lunar orbit, I began to 
review the theory presented in Christodoulidis et al. [1988], with frequent 
reference also to Kaula [1966]. In the midst of developing and de-bugglng my 
own programs for calculating these effects, communication with some of the 
Christodoulidis et al. co-authors (particularly S. Klosko and R. Williamson) 
led to their giving me a copy of the programs used In their work (it also led 
to the discovery, by me, of some errors in their programs...). By the start 
of the Spring 1991 AGU Meeting, I had also produced a list of ocean tide 
effects on the lunar orbit. These results, which are summarized in a 
manuscript currently in press as part of an IUGG Monograph (following the 
IUGG Meetings in Vienna), demonstrate that the bulk of tidal friction is 
produced by Just a few tidal constituents (M2, 01, N2, and also Mf and Ql). I 
found that my tidal theory failed to predict the observed amounts of Earth’s 
secular deceleration and changes in the lunar orbit, due to the theory’s 
abnormally small amplitudes for M2 and N2. However, the ability of any tide 
model to predict these effects in agreement with observation does not 
establish the general validity of that tide model — or conversely — since 
its prediction depends on so few constituents. 

At the Spring 1991 AGU Meeting, conversations with several people led me 
to reconsider all of the preceding work. Discussions in particular with D. 
McCarthy, R. Gross, and M. Eubanks emphasized that theoretical quantities such 
as the Earth’s angular velocity vector (and thus wobble and changes in the 
length of day) are not actually what is observed by current geodetic 
techniques, not even approximately so at high frequencies; for that matter, 
the very meaning of changes in the length of ‘day’ is unclear for excitations 
which possess semi-diurnal or diurnal periods. After the Meeting, as I 
attempted to understand this new viewpoint, which is especially important at 
high frequencies, I began to re-evaluate the Liouville equations themselves. 

As described in my latest manuscript, "Dynamic ocean tide effects on 
Earth’s rotation", which was recently accepted for publication by the 
Geophysical Journal International , the fundamental problems with the 
'traditionally formulated’ Liouville equations are Its failure to incorporate 
frequency-dependent fluid core decoupling (the correct formulation can be 
found in Smith & Dahlen [1981]!) and its failure to correctly model the 
response of the oceans to the excited rotational perturbations. In modifying 
the Liouville equations at this time, I took the opportunity to incorporate 
frequency-dependent Love numbers [Wahr 1981] and load Love numbers [Wahr & 
Sasao 1981] into both the Liouville theory and the spherical harmonic tidal 
theory. This in turn required that I carefully distinguish between prograde 
and retrograde phenomena, since the Love number frequency dependence stems 
from the retrograde nearly diurnal core resonance. 

Ultimately, then, I re-defined the tidal potential as 

U = Re{ U } U = M Y~ S e“ i<rt 

0 2 

for the tidal constituent of frequency <r, where M is the forcing amplitude 
and is the fully normalized complex spherical harmonic function of degree l 
and order n. This definition allows us to adhere to a rule that <r>0 
corresponds to prograde and <r<0 to retrograde phenomena. The generalization 
of my tidal theory to Include forcing by harmonics of negative order is easily 
accomplished: the only modification Is that in the matrixed momentum 
equations, eq. (6b) of Dickman [1989], 
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and 


are replaced with 


u 3 

a M vj 
o -* 



3 • it' and v * ll' , 

respectively, where it' is the collection of { } and 

jJi ( 0 for all (£, n) except (2, 

c * 1 -aH /g for 1 = 2, n = - s 

v o 


-s) 


The original Llouvllle equations are te.g., Lambeck 1980] 



where w = Clz + fiin Is the Earth's angular velocity vector and ^ Is the 
dimensionless excitation vector. The first equation allows us to determine 
the polar motion m s m^ + im^ produced by any Internal or external process; 
the second equation similarly yields the resulting axial perturbation m„ In 
rotation, thus the change In the length of day Alod * -86400nu sec and J the 
change in rotational angle Auti ■ - .fro^dt. Excitation Is possible from 

externally applied torques x, from internal changes t In angular momentum 
associated with currents relative to the rotating coordinate frame, or from 
perturbations c In the Earth’s inertia tensor. In these equations c. « A, C 
and « £1C; “that is, the perturbations In Inertia are small compared to the 
average principal moments of Inertia of the Earth (A = equatorial, C = polar), 
and the relative momentum perturbations are much less than the bulk angular 
momentum of the Earth. Finally, the equatorial components of the Llouvllle 
equation have been combined Into a single complex equation, with c * c +ic 
and 1 = !. + it , and the corresponding excitation function * ifi l + 1i/l 

Is multiplied by the scale factor <r r a Q(C-A)/A. 


The mantle, core, and oceans all respond to perturbations fin In Earth’s 
rotation. In order to determine the net effect of any rotational excitation, 
all of these responses must be characterized. In this report I will present 
their traditional formulations first and then, in stages, my modifications. I 
will go Into a fair amount of detail since the manuscript has not yet appeared 
in print. 


The mantle’ s response Is the most straightforward: if kE denotes the 
degree-2 Love number of the oceanless earth, then the changes In the Inertia 
tensor associated with the solid earth response are [Munk & MacDonald 1960] 


c 


kE 

ks 


(C-A)m 


'33 


4 kE 
3 ks 


(C-A)m 3 


where the ‘secular^ ^ove number (often called the ‘fluid Earth’ Love number) 
is ks » 3G(C-A)/(n a ), G is the gravitational constant and a is Earth's mean 
radius. (In the traditional view, the oceans' response is often included in 

kE. . . ) . 
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To the extent that the fluid core Is lnviscld, non-magnetlc, and axially 
symmetric, with polar moment of inertia C , it will be decoupled from axial 
perturbations in mantle rotation [Yoder et al. 1981, Wahr et al. 1981, 
Merriam 1982]; as measured from the rotating mantle, it therefore contributes 
a relative angular momentum 


At decade-scale periods, of course, geomagnetic torques will succeed in 
coupling core to mantle [Rochester I960]. 

Because the core boundary has the shape of an oblate spheroid, with 
flattening f , a wobble of the mantle engenders a slight response of the 
core fluid [Hough 1895]. As summarized in Wahr [1986], the relative angular 
momentum contributed by the core in this case is 


1 <r + (1+f A c“ +Q o' + (1+f )0 V ’ 

c c 

where A is the core’s equatorial moment of inertia (the first of these 
results may be obtained from Smith & Dahlen [1981] by combining their 
equations 6, 7, 10, and 11 in Appendix B with 4.11b). When <r is small 
compared to £1, this relation yields 1 proportional to m with no significant 
frequency dependence; however, in our work the frequency-dependent denominator 
in this expression cannot be treated as constant. 


Mantle elasticity allows the core-mantle boundary to deform slightly in 
response to the core flow [see Sasao et al. 1977, also Sasao & Wahr 1981]; the 
resulting products of inertia can be written 

* o' ~ i • 

c = <r + (i+f )n V = 0- + (i+f in V 

c c 


where £ = 2. 252x1 0" 4 and £s(A/A c )£. 

Combining the preceding mantle and core effects, the Llouville equations 
can be written 
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might be viewed as a frequency-dependent "Hough number" [cf. Dlckman 1983] and 
where now ♦ and denote additional excitations beyond those associated with 
the response of the oceanless earth to rotational perturbations. 
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Incidentally, the free wobble frequency <r for an Earth without oceans 

is easily found from these equations by setting = 0; since <r « Q, ic « 

(l+t)[l+f ] 1 « l+£-f and we have 6 

c c 


<r 
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<r r (l 


kE . 

ks 1 


* <V« ,A c 


(C-*) -g- 


As discussed in my manuscript, the oceanic effect on wobble period — 
determined using my most recent spherical harmonic tide model (including the 
features discussed presently) — implies that 2n/<r « 404.4 sidereal days, if 

the Chandler period is 434.4 sidereal days [Wilson 6 & Vicente 1980] ( — thus, 
an estimate of 30.0 sidereal days for the lengthening of the Chandler period 
by the oceans); this in turn leads to kE £ 0.31189. 

The free core nutation causes solid-earth Love numbers and load Love 
numbers to be strongly frequency-dependent at short periods [see, e.g., 
Melchior 1983], if the forcing potential is of harmonic degree 2 and order |s| 
= 1. Following Wahr [1981] and Lambeck [1988, Table 11.6] we first write 

<r - <r 

kE(<r) = 0.298 + — (1.23x10 ) (s = -1) 

<T “ <T 


where <r is the frequency of the 01 tide, chosen as reference, and <r Is 
the free core nutation frequency; <r = -2ir(l . 0021714) and from Lambeck 
[1988], <r Qi = -2rc(0. 926998) rad/sidereal day. Wahr [1981] also found 

kE = 0.299 (s * 0) 

kE = 0.302 (s ■ -2) ; 

these differ from the non-resonant portion of the s=-l kE as a result of 
mantle ellipticity and Earth’s prograde rotation. 

The oceanic response to polar motion and changes in l.o.d. will later 
require that we consider various types of forcing potentials; to keep track of 
the mantle’s response to all of those forcings we denote the Love number by 
kE , kE , and kE when the forcing potential is of order jsj =0, 1, and 2. 
Additionally, as 2 a result of mantle anelasticity, all Love number values must 
— as noted above — become equal to 0.312 at long periods, in order to 
produce the observed Chandler wobble period. Including all these 
considerations, we write the frequency-dependent potential Love number as 
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0.302 


s=-2 , 


kE (<r) = 

2 


0.312 


2n/|<r| s 4 days 
s=-2 , 2n/\<r\ > 4 days 


with 4 days chosen arbitrarily as the demarcation between elastic and 


and 


►onses 

of the mantle. 

We also define 

kE p (<r) 

= kE^cr) 

< r fc 0 


= kE^cr) 

<r < 0 


= k(<t) 

or fc 0 

% (<r) 

- k(<t) 

<r < 0 ; 


these definitions allow us to keep track of prograde versus retrograde 
forcing. 

Finally, quantifying the oceanic effects on rotation, as we will describe 
shortly, requires us to use various load Love numbers. As discussed in Wahr & 
Sasao [1981], the only load Love numbers possessing a core-resonance frequency 
dependence are those associated with the solid-earth response to loading by 
degree 2, order |s| =1 loads. From Wahr & Sasao, but noting also that Love 
numbers and load Love numbers should be equally dispersive, we write 


k'(<r) 
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-0.3083 


-0.3092 - 
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-0.3083 
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s=0 , 2n/|<r| a 4 days 

s=-l , 2w/|<r| a 4 days 

s=-2 , 2n/|<r| a 4 days 

all s , 2w/|<r| > 4 days 


where the long-period value for k' has the same scale factor, (0.312/0.298), 
as the long-period kE (and, as wfth kE, the scale factor is applied to the 
s=-l value of k' ) . The constant terms in these expressions have been altered 
slightly (by about 1/454) from Wahr & Sasao, to account for the slight 
differences between earth models 1066A and 1066B [see also Dahlen 1976] . 

It is useful to define k' = k' <<r) and k' = k' (<r) as the values 

202 202 2 1 2 \ 

of k' (<r) corresponding to loads of harmonic order s=0 or 2 and |s|=l, 
resplctlvely; and to distinguish 

k' (<r) = k' (<r) <r * 0 ; k'(<r) = k' (<r) <r < 0 . 

P 21 R 21 


Completing the description of Love number frequency dependence, it 
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follows from the same references as were used for kE and that we can write 



0.606 

<r - <r 

s=0 , 

2n/\<r\ 

s 

4 days 

hE(<r) = * 

0.603 + 01 (2. 46x1 0 -3 ) 

<r - <r 

FCN 

s=-l , 

2n/\<r\ 

3 

4 days 


0.609 

s=-2 , 

2n/ 1 <r | 

S 

4 days 


(0.603M0. 312/0. 298) 

all s , 

2n/ 1 <r | 

> 

4 days 


and 



-0.9994 

s=0 , 

b 

CM 

3 

4 days 

(<r) = - 

-1.0048 - - (2.88X10 -4 ) 

cr - <r 

FCN 

s=-l , 

2n/ 1 <r | 

S 

4 days 


-0.9994 

s=-2 , 

2ir/|<r| 

3 

4 days 


(-1.0048) (0.312/0.298) 

all s , 

2n/ 1 cr | 

> 

4 days 


Incorporation of frequency-dependent Love numbers into our rotational and 
tidal theories requires that these expressions must be substituted, where 
appropriate, into the Llouville equations (for the rotational theory), and 
Into the a factor and L matrix in equations (9) and (10) (for the tidal 
theory). 


Now for the oceans’ response. The complexities incurred because of the 
oceans can best be understood by considering the rotational perturbations in 
detail. The polar motion generated by a periodic excitation may be 
elliptical, including retrograde as well as prograde components, so for 
complete generality we write 


= Me** 

p 


M e 
R 




henceforth also allowing for the possibility of a complex forcing frequency, 
i.e. periodic motion (Re{<?}) plus damping (fm(£>). In terms of these 
variables, the change W in centrifugal potential accompanying wobble can be 
expressed as 


[W is presented in slightly different (but equivalent) form in Dickman 1988]. 

When an excitation generates a change In the length of day, the 
accompanying change W' in centrifugal potential also acts as a tidal 
potential; in this case the potential can be written 
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W' 



_2 2 
-£J a 



[see Lambeck 1988] where we set 


m„ ■ M'expdff-t) + conjugate 
O P “ 

The response of realistic (zonally asymmetric) oceans to W or W' will 
Include spherical harmonic components of all degrees and orders, exciting 
further changes in all components of rotation; with a dynamic response, these 
further excitations will derive from relative momentum as well as inertia 
perturbations. The initial tidal excitations also Include momentum and 
inertia contributions. From our spherical harmonic solution for tide height, 
the changes in the oceanic Inertia tensor can be calculated in a 
straightforward manner: 


- -T / TT v 4 ci*;<«-»{i;‘ eJlt - 1^1* 


-i<r*tj 


e - 


J 


'33 


2 /4iT 
3*5 


p a* Re 




— k' (<r))T° 
a 2 ~2 

i 


- V5(U 


a 

— k' )T° 
a o “o 

l 


:]} 


[cf. Dlckman 1989], where p is oceanic density and the tide height 

coefficients correspond to tKe tide height ‘truncated’ to zero over the 

continents [see Dickman 1988]. The present formula for c_« corrects a 

typographical error in Dlckman [1989 Appendix 2] mistakenly Involving 1/V7. 

In these formulae we have also Incorporated the change in mantle Inertia 

produced by ocean loading, via the load Love numbers k^, k^ [Lambeck 1980, 

1988], now treated as frequency-dependent. In the present work we also include 

a second consequence of core fluidity [Merriam 1980], which is that axial 

de-coupling of the core modifies the load Love numbers; in the notation of Nam 

& Dickman [1990], the modification factor Is a /a * 0.7919. 

3 1 


Defining 


r p (<r) 


v 4(i * k p ( ”' )) • * ■t'tF v‘ ( “ k ; (,)i • 


7 ' (o') 


■ *(i + __l k ' («r)) , 


and 


3 5 

a 


a 202 

l 


#co 2 (<r) * >^(1+ — )/(!+ — (<r)) t 
a o a 202 

l l 


we can write 


(17a) 

and 

(17b) 


, w -i I<rt , UT i,« -i<r*t 
c = r p (o , )I 2 e - - r R (<r) [J^] e - 


i<rt -i«r*t 

c e - + c e - 

p R 


C 33 = y 


' (<r)/?e|[T^ - K 02 (<r)T^]e'*-^j- . 


24 



From the solution for the tidal currents, the equatorial relative angular 
momentum of the oceanic response to W or W' can be calculated according to 


(18a) 

here 


. i<rt . , -l<r»t 
1 e - + 1 e - 

p R 


1 3 r* . A V' r 1 n nftk n 1 

2 P w a L 2^ L " ~ l Q *</ +1 > " ~ l C i Q iC/-l) J 


and the expression for 1 is the same as that for l p but with the tid£ velocity 
coefficients , v^ R replaced by their complex conjugates and Q replaced 

by (-l) n Q f “ n)4 . 1 C A , t ?, and Q™ are defined in Dickman [1988 and 

1989]. The formula for*l j^resented^ere differs from that listed in Dickman 
[1989]; the latter was actually the "motion" term (1 - 1c) of ¥ [Lambeck 
1980]. For consistency with the core and mantle effects, where 1 and c terms 
are explicitly treated in our Liouvllle equations, the angular momentum term 
here is restricted solely to 1. 

For the axial component of the relative angular momentum, we find 

(18b) l 3 = P^Keje^ 1 [ v J [ E* ]*} 

as in Dickman [1989]. Note that l involves harmonic coefficients of h rather 
than h. 


A traditional approach to the Liouville equations would Include the 
oceanic responses to W or W' through a standard Love-number formalism, for 
example writing c as directly proportional to m (or to n and to m*) with the 
proportionality constant an effective Love number. In the course of this 
project I was able to demonstrate, using my spherical harmonic tidal theory, 
that such an approach is not valid for a dynamic response of the oceans. 

After much effort, I discovered a relation between c and m which would remain 
valid under all conditions, including non-global oceans and other situations 
where the resulting polar motion may be elliptical and may be "cross-coupled" 
[Dahlen 1976] to the axial spin. This general relation can be written 


(19a) 

where 


c = Em + i Em + Fm + IFm + Gm^ + iGm^ 


E = 


m(<r)( B i l ] + <ry (<r) [ jT 1 ] 
- P - R BT2 

<ri + a 


E = 


R R^ P P~2 
<T» + O' 


F - - 


<r*r(<r)[ T*]* + <rr(<r) [ 

p fr^2 - R P^2 

<r* + <r 


F = 


7 (<r) [ T*]*~7 (<r ) l T^l* 
p m r n 

<r* + <r 


G = 


£%(«*)[ ( P i 'j; 1 ! 


^(cDKpI')*]* 


<r» + <r 


and 
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G 


-r p (<r)[( pi'); 1 ] - r^UT)^* 

<F* + <r * 

in these expressions the prograde and retrograde tidal responses to the 
potential W are given by 

f = M( J) and f = H(f) , 

P P n n 

whereas the tidal response to the potential W' is given by 

f = M' ( f' ) , 

p p 

and all tide height harmonic components in these formulae are ‘truncated’ . 

In a similar fashion, it is possible to write the relative angular 
momentum 1, and the axial perturbations to inertia and momentum c^» ly 
produced by W and/or W' , in terms of the rotational variables; the results are 

(19b) 1 = E n + JE b + F b* + IF ■* + G n„ + iC m„ 

L L L L L si L o 

(19c) c 33 = E'm + iE'm + F'ra* + iF'm* + G'm 3 + iG'm 3 

(19d) = E> + iE'm + F'm* + iF>* + G'm_ + iG'm_ . 

O L L Li L L O L O 

The various E, F, and G coefficients, which depend on <r, load Love numbers (in 
the case of c or c„,.), the parameters of the ocean model, and the tidal 
coefficients generated by W or H', are similar in form to those presented for 
c. 


The reader should note that determination of the E, F, and G coefficients 
requires a tidal theory from which solutions can be obtained to both tesseral 
and zonal forcing potentials (W and W' ) , at each tidal frequency. These 
coefficients cannot be estimated from known tidal solutions since they 
represent oceanic responses to potentials which are ‘structurally’ different 
than the original tidal forcing. And, tidal theories based on "hydrodynamic 
interpolation" are incapable of estimating the coefficients since there is no 
tidal data to which these responses can be tied. 


Finally, with periodic excitations, 
* 


. i<rt . -i<r*t 
¥ e - + * e - 

-p -r 


■ re 1 -* * (*')*e 


-P 


“P 


and equations (19) describing the oceanic responses to the Liouville 
equations can be reduced to 


( 20 ) 


D 11 M P + D 12 M r + D 13 H ; 


D 21 M p + D 22 M r + D 23 M ; 


-<r * 
r-p 


r ¥* 
r-R 


D 31 M p + D 32 M r + D 33 M ; 


-P 


These equations may be termed "broad-band" Liouville equations, because of 
their applicability at high as well as low frequencies. The D coefficients 
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are found according to 


A A 


'11 


22 


33 


'12 


21 


'13 


23 
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where 



fj_ kT p 

(<r). 

f, A A C x 

<r kT (<rk 

r p J 

( ? 
_ 0-1 + 

1 E L +£ * 

as — <T 

r | 

ks 

J 

+ ^ K P A 

a ks J 

-l 1 £1 

J A 


kT „ 

(cr)> 

* f, A A C 

sr 1 y* ... 

<r kT («•)•, 
+ r " 

I* + <rfi- . 

A. A 

- H V 1 * I 

= ~<T 

r 

l kS 

si 

Oil ^ , 

- [ R A 

a ks 

J 

j T V I I 

n Jl a J 

« 1 - 

c 

c + 

4 

kE (<r) 
C-A 0 

c '*i e l 

* 1 

— |T* 

G' 

L 

c 

3 

C ks 

C 

- C 
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E+-£-E 

kT (<r ) = kE (<r) + ks r — and kr (<r) 

p p A <r r R 


. , . . E+ a e l a 

kE (<r) + ks t — 

R A <r 


are, as a result of the oceanic responses, comp lex- valued as well as 
frequency-dependent whole-Earth Love numbers. 

Once the D coefficients have been calculated, tidal effects on rotation 
can be determined by using the tidal model and (17), (18) to construct the 
tidal excitation and then Inverting equations (20) for the unknown rotational 
perturbations H , H , and M' . That Is, when solved Jointly, the broad-band 
equations (20) yield R the polar motion m and axial perturbation m 3 produced by 
the Initial tidal excitation V or ^ 3 acting at frequency cr on an Earth whose 
mantle, core, and oceans all feed back to modify the polar motion or change In 
l.o. d. . 
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As indicated earlier, results have been obtained for all 32 of the tidal 
constituents modeled in Dickman [1991]. These results were also re-phrased, 
following the work by Gross [1991], in terms of the quantities directly 
observed in modern rotational geodesy: the so-called ‘nutation’ of the body 
axis (celestial ephemeris pole) relative to space; and the angular spin 
Increment (change in UT1). The results suggest that diurnal changes in UT1 
produced by ocean tides may be marginally detectable in modern high-accuracy, 
high-frequency space-geodetic data. Diurnal effects on polar motion may also 
be measurable. The detectability of semi-diurnal tidal changes in UT1 is less 
certain, because of the great disparity [cf. Brosche et al. 1989] in the 
inertia and relative momentum excitations predicted by semi-diurnal tide 
models. 
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